Transverse oscillating bubble enhanced laser-driven betatron X-ray radiation generation

Ultrafast high-brightness X-ray pulses have proven invaluable for a broad range of research. Such pulses are typically generated via synchrotron emission from relativistic electron bunches using large-scale facilities. Recently, significantly more compact X-ray sources based on laser-wakefield accelerated (LWFA) electron beams have been demonstrated. In particular, laser-driven sources, where the radiation is generated by transverse oscillations of electrons within the plasma accelerator structure (so-called betatron oscillations) can generate highly-brilliant ultrashort X-ray pulses using a comparably simple setup. Here, we experimentally demonstrate a method to markedly enhance the parameters of LWFA-driven betatron X-ray emission in a proof-of-principle experiment. We show a significant increase in the number of generated photons by specifically manipulating the amplitude of the betatron oscillations by using our novel Transverse Oscillating Bubble Enhanced Betatron Radiation scheme. We realize this through an orchestrated evolution of the temporal laser pulse shape and the accelerating plasma structure. This leads to controlled off-axis injection of electrons that perform large-amplitude collective transverse betatron oscillations, resulting in increased radiation emission. Our concept holds the promise for a method to optimize the X-ray parameters for specific applications, such as time-resolved investigations with spatial and temporal atomic resolution or advanced high-resolution imaging modalities, and the generation of X-ray beams with even higher peak and average brightness.

with the electron relativistic factor γ , the betatron period β , Planck's constant h and the speed of light c. In the bubble regime, the electron deflection parameter K is given by shaped gas jet (blue). The laser evolution during the first density peak leads to off-axis electron injection during the following density downramp. Subsequent large-amplitude betatron oscillations (yellow) cause emission of intense X-ray radiation (purple). The laser pulse is filtered out by a thin Al foil. The electron bunch is deflected and characterized using a dipole magnet spectrometer. The X-rays are measured using an absorption-filter based spectrometer. The inset shows gas density measurements for a distance of y = 1 mm (green) and y = 2 mm (black) above the gas nozzle using a backing pressure of 250 psi, the betatron experiments were performed at lower densities. www.nature.com/scientificreports/ with the betatron amplitude r β and the plasma wavenumber k p = n e e 2 /ǫ 0 mc 2 , where n e is the plasma density, ǫ 0 the dielectric constant, e the electron charge, m the electron mass and θ the maximum angle of the electron velocity to the wiggler axis. For typical LWFA parameters (where K ≫ 1 ) the number of emitted photons per oscillation period and electron scales as N ph ∼ K . Both, N ph and E crit can be enhanced and controlled with a larger betatron amplitude. N ph can also be increased with a higher electron beam charge and a greater number of betatron oscillations, which we experimentally demonstrate here.

Results
Overview. We compare the betatron emission from different plasma density profiles, including an approximately 7-mm-long double-peaked "M"-shaped profile (see inset Fig. 1) to flat-top density profiles with lengths of 4 and 6 mm using the same laser parameters. The M-jet has an effective plasma length of approximately 4.4 mm consisting of two 2.2 mm wide peaks that are interrupted by a density depression of 2.3 mm. The betatron source is driven by laser pulses with an energy of 2.6 J on target, a FWHM duration of 34 fs, focused to a 15 µ m FWHM spot using an off-axis parabolic mirror with 1 m focal length (see Fig. 1). The laser-plasma interaction in the gas target causes the generation of relativistic electron bunches and intense X-ray pulses. The laser pulse is blocked by a thin Al filter, while the electron beam and X-rays above 2.7 keV are largely transmitted (see Supplemental Material). The electrons are deflected from the beam path using a dipole magnet (10.2 cm, 7.6 kG), which allows the simultaneous measurement of the electron and X-ray spectra in a single shot. The electron beam spectrum and divergence are measured using a phosphor screen located at a distance of 1.95 m from the jet. The X-ray spectrum is diagnosed by the transmission through a filter array and an X-ray CCD (see Methods).
X-ray spectra. For each jet the X-ray flux was optimized by adjusting the plasma density and the jet position relative to the laser focal plane. The brightest X-ray pulses with the 4 and 6 mm flat-top jets were generated for a plasma density of 5 × 10 18 cm −3 and with the M-jet for 1 × 10 19 cm −3 . The observed spectral intensity for the "M"-shaped profile shows a clear enhancement in comparison to the other jets (see Fig. 2). Since the X-ray divergence exceeds our detector solid angle of 11 mrad , we estimate the total number of X-ray photons in a shot using the electron beam divergence. The total X-ray beam divergence is given by the convolution of the singleelectron emission and the maximum electron beam deflection angle inside the bubble. For a wiggler or betatron source, the X-ray divergence is dominated by the electron beam deflection 28 , which in the LWFA bubble regime is equivalent to the electron beam divergence after termination of the plasma, i.e. θ e,beam ≈ K/γ 27 . The X-ray divergence in this case is approximately given by θ X,beam ≈ θ e,beam ≈ K/γ . The electron beam generated by the 6 mm jet has a horizontal divergence of 10 mrad, whereas that from the M-jets are transversally clipped by the spectrometer dipole magnet with an aperture of 37 mrad (see Fig. 3). Using these values in the laser polarization dimension and the solid angle of the detector in the other dimension as conservative lower limits, we estimate an increase in the number of X-rays generated by the M-jet to be more than a factor of 6 compared to the 4 mm  Figure 2. Measured X-ray Spectra. (a) Spectral intensity extracted from the filter transmission for a flat-top 4 mm jet (red, dashed), a flat-top 6 mm jet (green, dotted) and the "M"-shaped jet for an interaction distance above the nozzle of y 1 = 1 mm (blue, solid) and y 2 = 2 mm (cyan, dash-dotted) for single shots. (b) A fiveshot average. The fit parameters for the five-shot average are given in www.nature.com/scientificreports/ and the 6 mm round jets (see Table 1). The sharp edges of the filter shadows indicate that the signal originates from a small source point and is not due to Bremsstrahlung background from scattering of the electron beam on the chamber walls. We can estimate the critical photon energy from the electron beam parameters using Eqs. (1) and (2) and the betatron period, which in the bubble regime is given by β = √ 2γ 2π/k p 27 . Using the electron beam divergence to estimate K, we obtain for the M-jet a maximum E crit = 186 keV ( θ e = 37 mrad, γ = 1100 , n e = 1 × 10 19 cm −3 ), for the 6 mm jet a maximum E crit = 5 keV ( θ e = 10 mrad, γ = 490 , n e = 5 × 10 18 cm −3 ) and for the 4 mm jet a maximum E crit = 12 keV ( θ e = 24 mrad, γ = 500 , n e = 5 × 10 18 cm −3 ). For the 4 mm jet, this estimate agrees reasonably well with the experimentally determined value. In case of the 6 mm jet, the spectrum extends to higher photon energies well beyond the estimated critical energy. This is likely due to the X-ray emission from electrons with higher energies at some point during the interaction, which occurs for a jet that is longer than the laser depletion length and the dephasing length (see section Discussion). The experimentally measured values for the critical energies of the M-jet are significantly lower, which is a result of the limited sensitivity of our spectrometer of 2.7-30 keV. Therefore, we expect a substantial fraction of the generated betatron photons to be outside our spectral detection limits, particularly in case of the M-jet. Considering the photons outside our angular and spectral acceptance, we speculate the actual enhancement in photon number to be even significantly larger for the M-jet to compared to the flat-top profiles.
The X-ray parameters can be modified through the plasma density profile. To this end we increase the distance between the laser-plasma interaction point and the gas nozzle. As can be seen from the density measurement in Fig. 1, the density gradient and the peak-to-valley contrast can be modified by changing the distance between the nozzle and the interaction point. We compare our measurements with a distance of y 1 = 1 mm to y 2 = 2 mm at similar plasma densities by compensating the gas jet backing pressure. The larger distance leads to a shift of the spectrum towards lower photon energies and an overall lower number of photons (see Fig. 2 and Table 1).
Electron spectra. We simultaneously measure the spectra of the electron bunches with angular resolution in the plane of the laser polarization (see Fig. 3). Although the spectra for all the jets have a broad energy distribution, their charge and angular distributions are significantly different. In particular, for the 6 mm flat-top density profile, we observe the smallest beam divergence, lowest charge and lowest peak energy (see Table 2). In stark contrast, for the 4 mm flat-top and the M-jet we observe a significantly larger divergence, which in case of the M-jet is even clipped by the angular acceptance of our dipole spectrometer of 37 mrad. From the electron beam divergence it becomes clear that only a fraction of the generated photons are observed by our X-ray camera that has an 11 mrad acceptance angle. While the 4 mm flat-top and M-jet both have large divergences, the total charge is significantly larger in case of the 4 mm jet. The maximum electron energy of the M-jet is significantly higher than that of the 4 mm jet, which is significantly higher than that of the 6 mm jet. Both, the 4 mm and the M-jet show energetically narrow features with large divergence. These features indicate off-axis electron injection at a specific position in the gas jet and the performance of coherent betatron oscillations throughout the subsequent density profile until termination of the gas jet. The large divergence, high charge and comparably low beam energy of the 4 mm jet in combination with the comparably weak X-ray emission indicates off-axis electron injection near the end of the gas jet during the density downramp. In case of the M-jet, the lower charge compared to the 4 mm jet, a higher maximum electron energy and stronger X-ray emission suggests off-axis injection during the density downramp of the first jet followed by subsequent acceleration and X-ray emission over a longer distance and possible loss of electrons due to the longitudinal bubble dynamics during the second density peak. Both, the lower maximum energy and lower charge of the 6 mm-long gas jet indicates Table 1. Gas jet comparison. The angular flux density is averaged over 5 shots for each jet. The total number of photons is estimated from the angular flux density and the beam divergence. As the X-ray beam divergence exceeds the detector solid angle, the horizontal X-ray beam divergence is inferred from the electron beam divergence and the vertical divergence from the angular acceptance of the CCD detector (11 mrad). The averaged critical energies E crit,i and amplitudes A i for each jet are obtained by fitting an N = 2 sum of synchrotron spectra (Eq. 3) to the 5-shot averaged spectra. The uncertainties are obtained by fitting a similar term to the 5-shot confidence bands.  www.nature.com/scientificreports/ electron injection and severe laser depletion well before the end of the gas jet, which leads to loss of electrons and decrease in beam energy due to dephasing. The stronger X-ray emission compared to the 4 mm jet despite a lower charge indicates earlier electron injection and betatron emission over a longer distance. These dynamics are qualitatively well in agreement with the dynamics that we observe in our our particle in-cell (PIC) simulations (see section Particle in Cell (PIC) Simulations). The electron and X-ray parameters of the M-jet can be modified by the density profile. Specifically, the steepness of the gradients in the modulated density profiles of the M-jets does not significantly impact the divergence and maximum beam energy, but the profile with less steep gradients decreases the accelerated charge, which has also been observed in numerical studies 29 . Due to these specific dynamics for each jet, the measured electron energies and charge after the termination of the jet have only limited correlation with the X-ray emission. Nevertheless, the measured divergence is indicative of the transverse electron momentum at the exit of the jet,   Table 2), while the divergence of the 6 mm jet (b) is more typical for LWFAs operated at unmatched plasma densities. The spectrum generated by the M-jet (c) extends to higher energies with a significantly larger divergence. The observed electron spectra, in particular the structured narrow-band, largedivergence features around 570 MeV agree well with the coherent betatron oscillations observed in our particle in-cell (PIC) simulations (see Fig. 4). The divergence of the electron beam is clipped by the 37 mrad angular acceptance of the dipole magnet. The spectrum of the M-jet using an increased interaction distance from the gas nozzle (d) has a similarly large divergence but less charge. For better visibility of the details, the plots are plotted on a logarithmic scale. A linear plot of the angular-integrated spectra is shown in (e). Table 2. Comparison of the generated electron beam spectra shown in Fig. 3. The third column is the charge integrated over an energy range, such that the lowest and highest generated critical photon energies lie in our detection range of 2.7-30 keV. to better understand the electron dynamics of the betatron source. In the simulations, the gas jet was represented by a combination of two Gaussian profiles that were fitted to the experimentally measured profile as shown in  The laser is moving to the right (i.e. in positive z-direction). The asymmetric longitudinal laser pulse profile leads to a transverse asymmetry in the bubble shape and off-axis electron injection during the density downramp (a)-(c). The field of the laser front is depleted by a half cycle, leading to a field spike with opposite sign and an asymmetry into the opposite direction (d)-(f). The controlled off-axis injection leads to a correlated transverse momentum distribution (c), (f), and (i) and coherent electron oscillations in the laser polarization (x, z) plane. Near the plasma density minimum, the bubble evolves into a super-structure, highly suppressing longitudinal accelerating fields (g)-(i). This leads to electron propagation with nearly constant energy as can be seen from the trajectories of 1429 macro particles from near the head of the electron bunch (k). The laser electric field E x is normalized to E 0 = mc 2 k p /e , the longitudinal and transverse electron momenta (color coded) are normalized to mc and the plasma density to the peak density n 0 = 1 × 10 19 cm −3 . www.nature.com/scientificreports/ the injected electrons perform collective large-amplitude betatron oscillations during the subsequent propagation (Fig. 4). In the laser conditioning section during the first density peak, the laser pulse undergoes substantial self-steepening and pulse compression, which leads to a nearly step-like temporal field distribution of the laser front ( Fig. 4a) 31 . This field distribution causes a transverse asymmetry in the bubble shape (Fig. 4b) 22 . The local depletion of the laser front 32 results in oscillating positive and negative amplitudes of the leading laser-field spike during further propagation. This leads to a bubble asymmetry in the laser polarization direction that changes direction in lock-step with the reversed sign of the leading field amplitude. As a result, the bubble boundary and in particular the large plasma density peak at the back of the bubble perform transverse oscillations in the laser polarization plane (Fig. 4). Due to the longitudinal bubble expansion during the density downramp in the injection section, electrons from this transversely oscillating high-density peak are injected into the accelerating fields at the back of the bubble. Because of the off-axis position, the electrons have a large correlated transverse momentum when injected, which leads to coherent large-amplitude transverse betatron oscillations in the laser polarization plane as they propagate (Fig. 4k) 22,23 . The longitudinally-correlated transverse momentum of the electrons can be seen in Fig. 4i. Despite a large transverse emittance integrated over the whole bunch, the local transverse slice emittance is small. Near the plasma density minimum, the combination of the laser evolution, increased plasma wavelength, lower plasma fields and injected charge (beam loading) leads to an elongated bubble super-structure, extending over multiple buckets (Fig. 4h and i). As a result, the longitudinal bubble fields are highly suppressed and the kinetic energy remains nearly constant for a large fraction of the oscillating electrons over an extended spatial region along the jet. The quasi-constant electron energy and coherent oscillations over a large distance can be seen in example trajectories of 1429 macro particles picked from a region near the front of the electron bunch (Fig. 4k). Due to the lower bubble fields during the density depression at the center of the M-jet, the electron beam experiences a transverse expansion and a small increase in the betatron oscillation amplitude around 3.5-5 mm 20,33 . However, this is not the main cause of the large-amplitude oscillations and enhanced betatron emission. During the subsequent density upramp of the second peak, the plasma wavelength decreases and a single-bucket bubble structure is slowly reestablished. As a result, electrons performing betatron oscillations are mainly confined to the first bucket, while electrons further back are lost (see movie in SM). Due to the increase in the bubble phase velocity, only a small fraction of electrons at the head slowly loose energy as they are propagating in a decelerating phase, while a large fraction of electrons are accelerated and gain energies similar to that of the bunch front. This leads to high-brightness X-ray betatron emission from electrons that perform coherent large-amplitude betatron oscillations with almost constant kinetic energy over a long propagation distance in the radiator section. Eventually, the laser pulse is nearly fully depleted and a bubble is driven by the injected electron bunch. The generated radiation for this density profile is significantly enhanced compared to flat-top gas jets. Both, the experimentally observed X-ray and electron beam properties agree well with the simulated kinematics of this process.

Scientific Reports
Jet comparison. We have also simulated the comparison between the different jet density profiles. The charge in an area near the laser propagation axis ( −15 < k p y < 15 ) and with a beam energy of γ > 20 is shown as a function of the laser propagation distance (Fig. 5). For the M-jet (blue), the copious amount of charge that is injected during the downramp after the first density peak has a maximum near the density depression and is subsequently slowly lost during the second density peak as the laser gets depleted. For the 6 mm flat-top jet (green), charge is injected after some propagation of the laser. The self-focusing and self-evolution during this propagation leads to an elongation of the bubble that results in further injection. During the subsequent propagation, charge is slowly lost due to laser depletion. In contrast, for the 4mm flat-top jet, charge is mainly injected during the density downramp at the exit of the jet. During the propagation through the jet, the laser self-evolves to a highly asymmetric longitudinal shape similar to that of the M-jet after the first density peak. Correspondingly, this leads to injection with a large transverse momentum.

Discussion
The measured X-ray and electron spectra indicate different electron beam dynamics in the three jets. The 4-mm flat-top jet generates electron beams with the highest charge and a large divergence while emitting the weakest on-axis X-ray flux. This indicates only a very limited distance during which radiation is generated, which suggests that the electrons are injected near the jet exit at the density downramp. This also agrees with an estimated acceleration distance of approximately 1 mm that is required to reach the observed electron energies for a plasma density of n e = 5 × 10 18 cm −3 . During the propagation through the jet, the laser significantly evolves and the large electron beam divergence indicates that the electrons are injected into a transverse asymmetric bubble driven by such a pulse. After injection, the electrons emit radiation during approximately only 2 betatron oscillations. The reasonable agreement of the X-ray spectrum with the estimated critical energy deduced from the measured electron bunch (see section X-ray Spectra) suggest that most of the betatron radiation is emitted from electrons with energies close to their measured spectrum at the exit of the gas jet and that the electrons do not undergo substantial dephasing. The electron beams of the 6-mm flat-top jet have the lowest measured charge and smallest divergence while generating a significantly higher X-ray flux than the 4-mm jet particularly at lower photon energies. The measured X-ray spectrum extends well beyond the critical energy of 5 keV, estimated from the measured electron beam parameters. This suggests that the betatron radiation is mostly emitted during the interaction from within the jet. In this case, electrons are continuously injected into the bubble 31 with a comparably large transverse momentum that is however significantly smaller than in the case of a self-evolved laser. The length of the jet is longer than the dephasing and laser depletion length, which results in a reduction of the electron energy and a loss of charge due to scattering out of the bubble. As the X-rays are generated over a longer distance throughout the electron propagation distance, this leads to the emission of X-ray pulses with a higher flux and higher photons energies than what is expected from the measured electron spectra at the jet exit.
The measured electron beams generated by the M-jet have the largest divergence with a charge that is less than that of the 4 mm-jet, while emitting the highest X-ray flux. The observed electron and X-ray spectra agree with the dynamics described in section Particle in Cell (PIC) Simulations, where the laser self-evolves during the first density peak into a highly asymmetric longitudinal shape. During the downramp of the first jet, electrons are transversally injected and X-rays are mostly generated throughout the density depression and the second density peak. Some of the injected charge is likely lost during the interaction. The M-jet combines the features of the other jets, namely the transverse injection of the 4 mm jet with the long emission length of the 6 mm jet for X-ray generation that is enhanced by the large betatron oscillation amplitude and by the density profile that allows X-ray emission over an extended length. While in the experiment it was not possible to measure the effect of only one of the density peaks, the 4 mm jet (at lower density) shows a qualitatively similar behavior to that of only the first density peak. These dynamics are consistent with simulations (see Figs. 4 and 5).
Increasing the distance between the gas nozzle exit and the laser interaction leads to a change in the plasma density profile, in particular increasing the width of the peaks, decreasing the density slopes and decreasing the peak-to-valley density ratio. In the experiment, this leads to a decrease in both, the injected charge and the on-axis X-ray flux. Simulations have shown that the slope of the density downramp has a strong effect on the electron beam parameters, in particular the amount of injected charge, the transverse emittance and the bunch duration 29 . This demonstrates that by varying the density profile, it is possible to modify the X-ray parameters.

Summary and conclusion
We experimentally demonstrate the novel Transverse Oscillating Bubble Enhanced Betatron Radiation (TOBER) generation regime that enables significant enhancement of laser-driven betatron radiation. In particular, we observe the generation of significantly more X-ray photons compared to standard targets using the same laser parameters. The X-ray source parameters are enhanced through manipulating the amplitude of the transverse betatron oscillations by using a laser pulse with a highly asymmetric temporal shape that has a nearly step-like rise. This leads to off-axis electron injection and coherent betatron oscillations, similar to those in a permanentmagnet wiggler. The source can be decomposed into three sections: (i) laser pulse shaping, (ii) electron injection and (iii) radiation generation. We demonstrate that this can be realized using a tailored plasma density profile. While longitudinally tailored density profiles have been used to improve the electron beam quality [34][35][36] , here we use a specifically tailored profile to significantly improve the X-ray beam quality. In our proof-of-principle experiments we use a profile that consists of two peaks separated by a density depression. Here, the laser propagation in the first peak leads to the asymmetric temporal pulse shape through self-evolution. The off-axis electron injection is initiated during the subsequent density downramp. The radiation is mainly generated during the density depression and in the second density peak. We demonstrate that the X-ray pulse properties can be modified by changing the plasma density profile.
The method can be extended to readily optimize the X-ray properties for specific applications. The plasma density profile enables control of multiple X-ray properties, including photon energy, number of generated photons and spectral intensity. The method also has the potential to adjust source properties, such as source size, beam divergence and likely pulse duration to the requirements of the specific application. Our concept has the potential for even further enhancement and control of the X-ray parameters and conversion efficiency. We expect that the X-ray beams can be scaled to higher peak brilliance through modifications in the plasma profile, by a higher-power driver laser or by using laser-generated electron bunches as drivers 37 to extend distance over which radiation is emitted. Conversely, the plasma profile can also be optimized to generate beams with higher average brilliance using lasers with lower peak power, operating at higher repetition rates.
The control and enhancement of X-ray parameters is of particular importance for the development of future variable and tunable compact, high-repetition rate X-ray sources that can be readily used in applications.

Methods
X-ray spectrum measurement and analysis. The X-ray spectrum is diagnosed by the transmission through a filter array consisting of Mylar, Al, V, Ti, Co, Fe, Zn, Cu, Zr, Mo, Pd, Ag, Sn and an Andor iKon-L SO X-ray CCD. The camera is shielded by a thin Al foil to minimize any optical background contamination. A Pb filter is used to determine any background. The filters were placed 1.9 m from the source and 0.75 m from the detector. The beamline filters and the quantum efficiency of the detector limited the spectral observation range to 2.7-30 keV (see Supplemental Material). The X-ray spectrum was deduced from least-square fitting of the transmitted signal through each filter, taking into account the transmission through additional beamline filters and the quantum efficiency of the detector (for detailed description see below). To account for the broad electron energy spread and the changing energy of the emitting electrons during the acceleration process, the spectral intensity was constrained to a sum of synchrotron spectra given by 28 where is the solid angle, E the photon energy and F is the modified Bessel function of the second kind. The amplitudes A i and critical energies E crit,i are fit parameters. We limit our model to two sub-spectra ( N = 2 ) to avoid over-fitting the data while showing a significant improvement over a single spectrum. A large set of Monte-Carlo (MC) simulations was used in a bootstrap method 38 to determine the confidence intervals via a statistical approach. We extract the X-ray spectra from the measured absorption filter transmission using a forward-propagation fitting method. To this end we first process the camera images using a despeckling algorithm to eliminate hot pixels, reduce the measured constant background given by the Pb filter transmission and perform a background gradient normalization using the direct beam in-between filters to account for beam nonuniformities. Experimental measurements for each absorption filter material are extracted by averaging the transmitted signal over the corresponding detector regions taken close to beam center. Standard deviations within these regions, prior to image processing, are used for measurement uncertainty estimates with appropriate scaling to account for image correction. The X-ray spectra are extracted via a Forward-Fitting Monte Carlo (MC) method. As input spectra, we assume a model constrained to a sum of synchrotron spectra given by Eq. (3). We use a sum of synchrotron spectra to account for the broad distribution and dynamics of electron energies, betatron amplitudes and periods, which evolve during the interaction. We found that a sum of two spectra give significantly improved the fit over a single spectrum. Due to a correlation of spectra with different critical energies, we limit the model to a sum of two spectra to avoid over-fitting the data. We iteratively randomly sample the multivariate parameter space (spectral amplitudes A i and critical energies E crit,i ) and forward-propagate the corresponding spectra through each filter, including the beamline and camera efficiency. We then compare to the experimental data through a least-square fit and determine the best fit parameters by minimizing the χ 2 value. For each shot we perform the Forward-Fitting MC procedure on 2000 independent data sets. Shots with a χ 2 per degree of freedom value of χ 2 /N > 0.2 were excluded from the analysis. We analyze this large parameter space sample to evaluate the goodness of the fit using the bootstrap method 38 and to determine asymmetric point-wise confidence bands. Specifically, for each photon energy value we find the 1-sigma distribution around our best fitting model, which is indicated as confidence bands in Fig. 2. From the single shot statistics, we calculate the average spectra and confidence bands using a weighted statistical approach. Shot-to-shot statistics are determined by point-wise averaging and appropriate uncertainty propagation considering the five shots with the highest signal intensity. The energies of the fits were limited to a range of 2.7-30 keV.

Data availibility
The datasets generated and/or analysed during the current study are not publicly available due to size but are available from the corresponding author on reasonable request.